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ABSTRACT 


This thesis describes the results of an investigation of the throttl- 
ing effects of a transverse jet when injected into a sonic throat. The 
process, commonly referred to as aerodynamic throttling, was studied both 
analytically and experimentally. 

An analytical model is developed utilizing momentum considerations. 
This model successfully predicts the measured throttling results for 90° 
injection and mainstream flow reductions fom approximately 60 percent. 
The agreement between theory and experiment becomes less satisfactory as 
the upstream angle of injection is increased and at extremely high 
throttling conditions. 

An experimental program is described in which a two-dimensional 
nitrogen jet is injected into the mainstream through slots of various 
widths and inclinations. Visual flow patterns are shown to qualitatively 


describe the nature of the interaction between the jet and the mainstream. 
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NOMENCLATURE 


a Speed of Sound 


B Mainstream duct height 

Cy Discharge coefficient 

d, Width of jet nozzle throat 

ey Effective pressure force on "windward" side of jet; acting in the 
x-direction 

h Maximum height of sonic line above the- wal] 


H Jet slot width nondimensionalized by the mainstream duct height d,/B 
m Mass Flow Rate 

M Mach number 

E Pressure 

i Temperature 

V Exit velocity of jet 


x,y Coordinates along, or normal to wall, respectively 


Oo Angle of jet injection with respect to the vertical 
Y Ratio of specific heats 
9 Density 
SUBSCRIPTS 
j Jet conditions 
m Mainstream conditions 
G Reference conditions based on mainstream flow with no jet injection 
0 Upstream conditions for mainstream 
1 Mainstream conditions at point of injection 
Z Mainstream conditions at point of maximum jet penetration 





P(M) 


SUPERSCRIPTS 


Stagnation conditions 
Sonic conditions 
FUNCTION DEFINITIONS 
Area ratio in locally SIERO AICINTAOW 
] : = a BAIT 
= a 


Pressure ratio in locally isentropic flow 


= a 
I gly [7 
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I. INTRODUCTION 


Seventeen years ago an article was published by Manoury, et.al. 
[Ref. 1] proposing an exit plane aerodynamic nozzle for a jet engine. 
This article pioneered the general concept of aerodynamic throttling. 
Aerodynamic throttling is the process by which the mass flow rate of a 
primary flow is controlled (throttled) by injecting a secondary (jet) 
flow into the mainstream. The jet/mainstream interaction creates what 
is referred to as an “aerodynamic throat" and has the effect of reducing 
the mainstream mass flow rate. Since Manoury's article in 1955, the 
feasibility of aerodynamic throttling has been shown experimentally 
[Ref. 2-4] and predicted analytically [Ref. 5-7]. The knowledge gained 
from these studies has been used by designers of systems that involve 
the injection of one fluid into another such as variable-throat nozzles, 
combustion chambers, and rocket nozzles cooled by injection. 

This investigation was undertaken with the idea of applying the 
principle of aerodynamic throttling to a fluidic device which could be 
used in a missle control system. It was envisioned that a modified pro- 
portional amplifier, with throttling control, could be used as a moment 
producer for a missile steering system employing external jets. The 
advantage of such a device would lie in the fact that during the no- 
control phase of operation, the gas expended by the fluidic control 
System could be throttled thereby saving gas and allowing for either 
more control time, or a reduction in required control gas with a 
resulting increase in missile payload. 

An important aspect of the aerodynamic throttling process is the 


extent to which the jet to mainstream total pressure ratio affects the 
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penetration of the jet into the mainstream. This study attempts to improve 
the analytical prediction of the penetration height of the jet and also 

to extend the information available about the throttling process by 
experimentally determining the effect of high total pressure ratios. 
Previous investigators have not gone above jet to mainstream total 


‘pressure ratios of 5:1. 
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II. ANALYSIS 


A. GENERAL CONSIDERATIONS 

The goal of this study is to examine the jet penetration and throttl- 
ing characteristics of a transverse jet injected into a sonic throat. 
The jet, shown in Fig. 1, is sonic, underexpended, two-dimensional and 


inclined at an angle a into the mainstream. 
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Figure ]. Region of Interest. 

The basic structure of the jet in the vicinity of the nozzle exit is 
similar to that of a free jet expanding into a quiescent atmosphere 
except that the pressure in region A (windward side) is higher than that 
in region B (bubble side). This pressure imbalance results in a jet 


plume that is skewed relative to the nozzle axis. For a given injector 


ie 





geometry and specific gas properties, the jet total pressure and the 
exterior pressure in regions A and B determine the location of the jet 
boundaries. The jet flow in the core of the plume quickly becomes over- 
expanded and must pass through a strong jet shock while the flow along 
the boundary of the jet passes through a series of oblique shocks and 
remains supersonic. The strong jet shock, known as the Mach disk, is 
nearly normal to the jet flow direction. The jet flow immediately down- 
stream of the Mach disk is subsonic and at a pressure that is higher 
than the mainstream pressure. Beyond this point the jet reaccelerates 
and expands in an attempt to achieve pressure equilibrium with the 
mainstream. Mixing between the jet and the mainstream is considered to 


be negligible in the region prior to the Mach disk. 


B. JET PENETRATION 

In the jet penetration studies conducted by Schetz and Billig [ref 8] 
and Orth and Funk [Ref. 9] the characteristic dimension of the inter- 
action region was taken to be the height to the center of the Mach disk. 
This dimension was considered to be a measure of the jet penetration into 
the mainstream. Barnes, et.al. [Ref. 10] as well as Spaid and Zukoski 
[Ref. 11] chose as the characteristic dimension the height of the sonic 
line occurring in the jet at a short distance beyond the Mach disk. The 
actual penetration of the jet can only be defined in terms of its effect 
upon the crossing flow. The characteristic height, whether that of the 
Mach disk or the sonic line, is useful mainly as a scale for the pene- 
tration. Both of these heights have been found to be valid parameters 
for the correlation of jet interaction data. In this analysis the sonic 
line is used as a more convenient characteristic height and, in addition, 


it is felt that this measure is more nearly the actual or "effective" 
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depth of penetration as discussed below. To the extent that the calcula- 
ted penetration is the actual penetration, the throttling characteristics 
of the system can be calculated by a direct application of the principle 
of mass flow continuity. 

As the jet is turned by the bounded mainstream, the Mach disk becomes 
more nearly perpendicular to the wal]. In this configuration, a measure 
of jet penetration based on the distance to the center of the Mach disk 
would only account for the penetration of a portion of the jet. Simply 
adding half the width of the Mach disk to this height would still not 
account for the full penetration. Vick, et.al. [Ref. 12] has shown that 
for three-dimensional free jets at low total pressure ratios, the Mach 
disk occurs close to the point at which the jet plume achieves its 
maximum width and that the Mach disk diameter is no more than sixty 
percent of the maximum jet diameter. It appears reasonable then to 
assume that for a two-dimensional jet the maximum penetration into a 
bounded mainstream would be some distance greater than the height to 
the center of the Mach disk. 

Although the flow immediately behind the Mach disk is subsonic and 
the mainstream would be expected to stronly affect the flow in this 
region, there is jet flow along the windward boundary of the jet that 
has not passed through the Mach disk and remains supersonic for some 
distance downstream of the Mach disk. Vaughan [Ref. 13] and Maurer 
[Ref. 14] have observed, using shadowgraphs and oi] films respectively, 
that the jet achieves a sonic velocity within a very short distance after 
the Mach disk and that negligible mixing occurs between the jet and 
mainstream prior to the sonic line. Thus, it is hypothesized that the 


Supersonic flow along the boundary shields the subsonic core of the jet 


is 





from massive mixing with the mainstream until the region of the sonic 
line is reached. Beyond this point it is assumed that the jet loses its 
identity, mixes with the mainstream, and ceases to affect the mainstream 


in the sense of penetration or throttling. 


C, METHODS OF ANALYSIS 

1. Previous Work 

The problem of a secondary jet issuing into a bounded mainstream 

for the purpose of throttling the primary flow has been analyzed by 
numerous investigators in the past. In a recent analysis, Nunn [Ref. 15] 
makes use of a momentum balance across a control volume which includes 
both the mainstream and the jet. This momentum balance is then combined 
with continuity expressions for the mainstream and the jet along with a 
Statement of pressure equilibrium between the jet and mainstream at the 
point of maximum penetration. To account for jet total pressure losses, 
it is noted that the strongest shock that can occur in the jet is ata 
jet Mach number which corresponds to an isentropic jet expansion to the 
Critical mainstream static pressure. The solution of the system of 
equations presented in Nunn's work is based on the assumption that the 
Strongest shock (i.e., the Mach disk) does occur in the jet and that it 
is located at the point of maximum penetration. It is felt that this 
assumption places an unrealistic constraint upon the flow field when a 
wide range of total pressure ratios is considered. 

Nunn's analysis determines the shock strength purely on the basis of 


the ratio of total pressures: 
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This is equivalent to assumina that the jet is exhausting into a quiescent 
atmosphere ce = Pa and does not take into consideration the crossing 
flow and the effect of the turning of the jet. Due to the curvature of 
the jet plume, the flow within the jet experiences an additional compress- 
ion that will cause the shock to occur at a lower Mach number than would 
be predicted for a free jet. 

It has been shown by Adamson and Nicholls [Ref. 16] as well as Vick, 
et.al. [Ref. 12] that for free jets, the distance to the Mach disk follows 


an approximately linear variation with: 


O 
PE 
a 
P 


co 


when plotted on log coordinates. Additionally, Vick has shown that the 
ambient pressure gee has a definite effect on the Mach disk location. 


For a given pressure ratio: 

O 
J 
P 


co 


~~ 


the lower the ambient pressure, the less the distance from the nozzle 
exit to the Mach disk. Based on the above findings, it would seem un- 


realistic to assume that the Mach disk, for all pressure ratios: 


oO 
O |e. © 
w 


would position itself at the point of maximum penetration. 

It is advantageous to seek an analytical model whereby the jet bound- 
aries as well as the Mach disk are more accurately predicted for jets 
which experience not only a different pressure on each side, but also 
encounter a varying pressure on the windward face. A computer program 


was written employing the method of characteristics to determine the 
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shape of the jet plume along with the Mach number distribution on the 
centerline. This program was successfully completed for jets which 
experience a constant pressure on each face, but attempts to modify the 
program to account for a varying pressure were discontinued for two main 
reasons. First, it is not known a priori just what the pressure distribu- 
tions are on the faces of the jet, and secondly the complexity of the 
computer program that would be necessary to handle this problem by means 
of the method of characteristics was deemed to be prohibitive. There- 
fore, it was decided that an alternate approach to the problem was 
needed. For purposes of completeness, however, a brief description of 
the method of characteristics calculations, together with a computer 
listing, is included as Appendix A to this report. A primary considera- 
tion of this alternate approach was to avoid the assumption that the 
Mach disk occurs at any‘particular location in the jet. 
2. Present Analysis 

As indicated previously, the jet is assumed to be sonic at 
injection, underexpanded, two-dimensional and inclined at an angle a 
into the mainstream. The jet undergoes an initial isentropic expansion 
and is then recompressed by a strong normal shock, the Mach disk. The 
jet fluid leaving the Mach disk again expands to sonic velocities a 
short distance downstream. The flow is assumed to reattach itself to 
the wall and continue to flow parallel to it. 

The dotted line in Fig. 2 is designated as the sonic line and 
is the transition from sonic to supersonic velocities. The actual 
geometry of the sonic line is unknown. The geometry has, however, no 
effect on the results of the analysis if it is assumed that the turning 
of the jet occurs prior to the sonic line so that the average momentum 


of the jet fluid at the sonic line is parallel to the wall. This 
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assumption is the same as that made by Barnes, et.al. [Ref. 10] and 
Cassel, et.al. [Ref. 17] in their studies of jets in an unbounded 


mainstream. 





Figure 2. Schematic Diagram of Throttling Region Showing 
Control Volume. 


A control volume (dashed line in Fig. 2) is defined which is 
bounded by the windward face of the jet plume from the nozzle to the jet 
sonic line (this is also the point of minimum mainstream flow area), the 
jet sonic line to the wall, and the plane of the wall. 

Conservation of jet momentum requires that the change in the 
X-component of momentum be equal to the forces in the x-direction acting 
upon the control volume. Recalling that the momentum vector of the jet 
fluid at the sonic line is assumed to be parallel to the wall, the 


bd 


conservation of momentum can be expressed as follows: 
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F, ~ PS hem, [a + V, sin a] (1) 

In this expression Pe is the average pressure force acting in 
the x-direction on the windward face of the jet and will be evaluated in 
subsequent paragraphs. It has been further assumed here that the shear 
stresses acting along the jet/mainstream interface do not contribute 
Significantly to the x-component of the net force on the plume. This 
assumption is characteristic of earlier jet interaction analyses (includ- 
ing those previously mentioned) and is based upon the hypothesis that 
the drag due to Shear will only play a significant role at downstream 
regions where turbulent mixing is the dominant exchange mechanism for 
mainstream and jet momentum. This assumption presupposes an initial 
transverse jet momentum sufficient to lead to a large penetration in a 
Short distance. As the distance to jet maximum penetration increases 
(as, for instance the upstream angle of injection increases) the neglect 
of shear stresses becomes less credible. | 

Equating the mass flow rate at the nozzle throat to that at the 


sonic line yields: 





~ 2 2 % 2 =r Jot 
ms = p4 dy, a, = p hav = y(4y] (2) 


where the subscript (t) refers to conditions at the injection sonic 
throat. The latter relationship in Equation (2) restricts the analysis 
to conditions in which the jet behaves as an ideal gas and the process 
is adiabatic. Thus, for situations in which the jet total temperature 
varies widely from that of the mainstream, it may become necessary to 
consider the transfer of heat across the jet mainstream boundary. For 
the adiabatic process, a = a, and Equation (2) produces the following 


relationships: 
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It should be noted that Equations (2) and (3) are the same as those 
obtained by Barnes, et.al. [Ref. 10] in their study of transverse jets 
interacting with an unbounded free stream. Substituting Equations (2) 
and (3) into the momentum equation yields an expression for the maximum 


height of the sonic line above the wall: 


0 a. 
. >» \Y~1 Ve 
aa j 2 | ( : | 
aa <a y il +— sina]? I : (4) 
d, ie (7 ay 
h 


Since the experiments that were conducted for this study treated only 





jets which were sonic at injection, V, = a,, Equation (4) takes the 


e 
final form of: 


=] 
(2) [+ (1 + sin 9) +1], (5) 


The height of the sonic line, h, as given by Equation (5), is a character- 
istic dimension of the flow field which according to the previous 
arguments 1s an estimate of the depth of jet penetration into the 
mainstream. 

If Equation (5) could be solved directly for h, then it would be 
a simple matter to determine the effective throttling of the mainstream. 
The problem, however, lies in determining the effective average pressure, 
Fh, acting on the windward boundary. Since the mainstream is bounded 
on all sides, a sonic throat will be formed in the mainstream at the 


point of maximum jet penetration. The pressure acting on the windward 


bd 
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face of the jet is assumed to vary from the mainstream static pressure 
* 
at station 1 [Ref. 15] down to - at station 2. See Fig. 2. 
The effective average pressure acting on the windward boundary 


can be expressed in integral form as: 


2 
ee at (6) 


h yy 


Expressing this integral in terms of Mach number gives: 


Fe = Pao fo vim & ay (7) 
= hy dM 
mo 
Relating y to M 
Eye en et) 
B A, ACM) 


and differentiating y with respect to M gives: 


foi OM 


Upon evaluating & LA(M)] Equation (8) becomes: 


yt] yy 
dg . BP ll me 2 
dM AMY dM 1+ tl 


and substitution of this expression into Equation (7) gives the following 


for the effective pressure: 


0 yt 
Z- a - Svan AM) Jt 2d aw (9) 
Mn AM) JM 4 eh? 
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Finally, if Equation (9) is substituted into Equation (5) an equation in 


terms of the mainstream Mach number, thee appears as follows: 


+] 
sty +1 
eT | 2y-TY 
. 2 mo ieee x1 2] / | (2, + Se) Oe ees 
mo y+ M 2 ue y+ 2 
20 mo 
dy —L.- 
= i = y (1 + sin a) + 1 (2) (10) 
7 yt] 


0 


For a given set of parameters a a, and so Equation (10) can be solved 


for the mainstream Mach number Ms uy 


From the isentropic area relationship: 


> Bn 7 AMM) 


an expression for the maximum jet penetration, h, can be obtained: 


h A(M_) - | 
B | -AW,o) a 


D. THROTTLING CHARACTERISTICS 


A reference mass flow rate is defined as that which would flow 


isentropically in the mainstream with no jet injection; 


where 


2] 





Using this reference mass flow rate, the throttled mainstream and jet 


mass flow rates may be expressed in nondimensional terms as: 





m 

m | 

= = A(M.__) (12) 
m. A He 

m. d, p? 

— ~ (13) 
m. B _ 


The theoretical curves for flow throttling are shown in Figure 3. 
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Figure 3. Throttling for Various Injection Angles 
(Unmodified Theory). 
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ITI. APPARATUS AND EXPERIMENTAL METHODS 


The goals of the experimental portion of this study have been to 
evaluate the analytical model and examine the throttling effects of a 
transverse jet with various slot widths, dis and angles of injection, a, 
when injected into a sonic throat. To achieve these goals, it was 
desired to experimentally determine the mass flow rates of the mainstream 
and jet during the throttling process. 

The experimental apparatus consisted of four distinct functional 
units: the first being the primary air flow system; the second being 
the test section itself; the third is the jet or secondary flow system; 
and, the fourth is the instrumentation for gathering flow information. 

The primary system, shown in Fig. 4, consisted of a 200 psig air 
compressor and two 117 cu. ft. tanks. With both tanks pressurized to 
maximum capacity, it was possible to conduct blow-down tests of approxi- 
mately one minute duration with a reference mass flow rate of 1.19 Ibm/ 
sec. The air left the tanks and passed through a pressure regulating 
valve just downstream of each tank. The pressure regulators were 
remotely controlled by a compressed nitrogen control system and were 
capable of regulating the mainstream pressure from 0 psig to 150 psig. 
After the regulators, the flow from both tanks was combined in a single 
3-inch pipe in which was installed a remote controlled on-off valve. 
Downstream of the on-off valve, the 3-inch pipe was fitted with a 
Standard ASME flat plate orifice which was used to obtain mass flow 
measurements for the primary system. The air then flowed into a 5-inch 
pipe and finally to a transition section that changed the circular 


cross-section to a rectangular cross-section measuring 2.625" x 3.0". 
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Primary Flow System with Control Panel. 


Figure 4. 
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The 5-inch pipe and subsequent transition section were part of a small 
wind tunnel which had been obtained as surplus government equipment. 
Since the available compressed air tanks could not provide sufficient 
blow-down times through the existing cross-section, the test section 
that was installed in the tunnel had plexiglass inserts on either side 
which reduced the width of the tunnel from 2.625" to 1.4". In this way 
the rectangular cross-section just prior to the test section had the 
dimensions of 1.4" x 3.0". 

The test section, shown in Fig. 5, consisted of a countoured upper 
block and a flat lower block which contained the jet slot. The con- 
toured upper block choked the flow at a 1.4" x 0.737" throat. The air 
then exited through a nozzle which was divergent on one side only, the 
upper block; the straight wall of the lower block served as a plane of 
symmetry for simulating symmetric injection. The lower block was fitted 
with variable inserts to allow for a change in both the slot width, d,. 
and the angle of injection a (see Fig. 6). Care was taken in the design 
of the inserts to insure that the jet occupied the entire jet slot 
(from side wall to side wall) upon entry into the mainstream flow. 

The secondary flow followed a path similar to the primary system in 
reaching the test section. Eight bottles of compressed nitrogen, mani- 
folded together, provided the secondary flow. (See Fig. 7). After 
leaving the manifold, the nitrogen passed through a remote controlled 
pressure regulating valve. The regulator would reduce the 1800 psig 
bottle pressure of the nitrogen to values ranging from 0 psig to 1500 
psig. The nitrogen then flowed through a remote controlled on-off valve 


and on through a sectionof 1" stainless steel tubing which was fitted 
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Figure 5. Test Section 





Figure 6. Lower Block with Removeable Insert. 
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with a standard ASME flat plate orifice to a plenum chamber located in 
the lower block of the test section. From the plenum chamber the second- 
ary flow was exhausted through a slot of variable width, and into the 
primary flow at the sonic throat. 

The instrumentation for both the primary and secondary systems con- 
sisted of variable reluctance pressure transducers whose output was 
continuously recorded on strip chart recorders. In this way, the orifice 
pressure drops and corresponding upstream nemnee required for deter- 
mining mass flow rates by the ASME Power Test Code methods, [Ref. 18], 
as well as the total pressure measurements for both the primary and 
secondary system, were obtained as a continuous record: during any 
experimental run. The total pressure probe for the primary system was 
mounted in the transition section just prior to the test section, and 
the total pressure probe for the secondary system was mounted in the 
plenum chamber in the lower block of the test section. Along with the 
output from the pressure transducers, bourden tube gages installed in a 
control panel, Fig 4, registered the total pressure of both systems, 
reservoir pressure of the primary system, and control pressure for the 
primary and secondary pressure regulating valves. The remote controls 
for the reducing valves as well as the on-off valves were also located 
on this control panel for convenience of operation. 

The mainstream supply pressure, set at ia = 35.0 psig, was chosen 
for convenience and to provide sufficient run times. The flow of the 
unthrottled nozzle was measured for mainstream total pressures of 30.0 
psig to 50.0 psig with an average discharge coefficient in the range of 
1.05 + .02. The discharge coefficient used here is defined to be the 


actual flow divided by the theoretically predicted flow. Considerable 
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effort was expended to reduce the system leakages to a minimum value so 
that all of the measured primary flow actually passed through the throat. 
The value of the discharge coefficient given above indicates that this 
effort was not completely successful but further reductions would have 
entailed massive modifications to the equipment on hand. Further 
comments concerning the difficulties with the experimental apparatus are 
given in the section on recommendations for future work. At the supply 
pressure of 35.0 psig and air temperature of 530°R the resulting mass 
flow rate of the mainstream air was 1.19 Ibm/sec. This value was used 
for m, in the reduction of data after correcting for slight variations 
in eae 

The test section was provided with windows so that visual observations 
of the jet/mainstream interaction could be made. Visual flow patterns 
were obtained by ee inside surfaces of the plexiglass inserts 
with a mixture of titanium dioxide (Ti0.) and 200 Dow Corning Silicone 
Oil with a viscosity of 350 centistokes. The oi] and titanium dioxide 
were mixed in the proportions of 7cc of oil with Icc of T10.. The 
mixture was applied with a brush with little regard to the coating 
thickness since the viscosity is the parameter which dictates the speed 
which which the coating is blown off. Photographs were taken of the 
resulting flow patterns. Although it is difficult to obtain any 
quantitative results from these pictures, they were useful in 
qualitatively examining the jet/mainstream interaction region. 

A simplified linear error analysis was performed on the experimental 
data in an attempt to arrive at bounds on the uncertainty of the mass 


flow ratios of: 


= 


ls 


: e 


and 


=j e 
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The sources of error in the experimentally determined mass flow rates 
fall into three main categories: the first being the error associated 
with the various parameters which enter into the equation for: determining 
the mass flow rates from the ASME Power Test Code; the second is the 
error in the pressure transducers and associated electronics; and, the 
third is the error in reading the pressure traces on the strip chart 
recorders. 

The error associated with the various constants appearing in the 
ASME mass flow rate equation was found to be insignificant. Furthermore, 


when examining the error in the ratio of mass flow rates, i.e., 


it was found that any error in the previously mentioned quantities 
would have no effect on the error of the ratio since the same error 
appears in M and m, and would therefore cancel out. The total error 
in the ratio: ; 

um 


r 


was found to be a function of the error in measuring the pressure drop 
across the orifice plate, as well as in measuring the upstream static 
pressure. Furthermore, the error in both of these pressure measurements 
1s dependent upon transducer error and reading error. 

Each individual pressure transducer along with the associated 
electronics and strip chart recorder was calibrated as a total unit. In 
this way, much of the error which might have been inherent in single 
components was calibrated out of the system. It is felt that the 


transducers and electronics used for determining the upstream static 
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pressure and the total pressure for both the primary and secondary 
systems, were capable of determining the pressure to within + 1.0 psi. 
Combining the transducer error with the error associated in reading the 
trace on the strip chart recorders, provided a measure of the total error 
in determining the upstream static pressure and total pressure of the flow. 
It was found that the static pressure and the total pressure of the 
primary and secondary systems could be determined to within + 1.25 psi 
and + 3.0 psi respectively. The sensitivity of the differential trans- 
ducers used for “eastiing the pressure drop across the orifice plates 
was such that the differential pressures could be determined to within 
+ 0.1 psi. This value includes the error in reading the strip chart 
recorder. By using the previously determined errors for the individual 
pressure measurements, upper and lower bounds were determined for the 


uncertainty in the nondimensionalized ratios of: 


m m- 

= ance 

m. m. 
for: 

m 
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the percent uncertainty is + 0.3 percent, and for: 


ll 


O12 


m 


the percent uncertainty is + 24.0 percent. The high error for the case 


of: 


0.2 


(high throttling) comes from the measurement of the primary orifice 
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pressure drop. As previously mentioned, the total uncertainty in 

measuring the pressure drop across the orifice plate was + 0.1 psi. For 
low values of: ; 
“mn 


My 


(high throttling) the differential pressure is less than 0.5 psi, there- 
fore the percent uncertainty in this measurement alone is on the order of 
+ 20.0 percent. This problem could be corrected if the flow measuring 
device would give significantly larger pressure drops at the lower values 


of: 

Mn 

My 
In the case of an orifice meter, this would entail using orifice plates 
with smaller diameters as the throttling was increased. The experimental 
setup for this study did not lend itself to changing the orifice plates 
without Benondenabie effort, so all experimental runs were conducted with 
the same size orifice plate. It should be noted however, that the 
consistency of the data seems to indicate that the data falls well within 
the bounds of uncertainty. A similar analysis was performed for the 
ratio: 

Mm, 
and it was found that the percent uncertainty in this ratio was on the 
order of + 1.5 percent over the entire range of: 

=a 


m 
r 
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Whereas the uncertainty is given for the ratio of mass flow rates, 
it must be remembered that the uncertainty for the primary or secondary 


mass flow rate alone, may differ significantly from the uncertainty for 
the nondimensionalized ratio. For example, if: 


m 
= = 0.484 


m 


the percent uncertainty in this ratio iS + 3.8 percent. 


However, the 
percent uncertainty in m alone is +5.5 percent. 
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IV. DISCUSSION OF RESULTS 


Throttling measurements were obtained for values of H of 0.02, 0.04 
and 0.08 for perpendicular injection (a = 0°). Measurements were also 
taken for a = 15° with H = 0.023, and a = 30° with H = 0.027. For each 
of these five settings the jet supply pressure was varied so that data 
were obtained for: 

il 

my, 
over a range of 0 to 0.6. Figures 8 and 9 show the results of the 
throttling measurements. 

Examination of the data in Fig. 8 shows that although the data is 
consistent and repeatable, there is a general departure of the theory 
from the data. An improvement on the analytical model will be discussed 
in detail below. 

In Fig. 9, it 1s evident that the effect of angle of injection as 
predicted by the theory is qualitatively verified. However, it appears 
that as the angle of injection exceeds approximately 30° the effect on 
the throttling becomes negligible. This phenomenon has also been observed 
by Nunn [Ref. 6]. At high angles of injection it is probable that the 
primary air flow partially fills the jet slot and the effective angle of 


injection becomes something less than the designed-for value. (See Fig. 10) 


The departure of the theory from the data at high angles of injection 
can be partially attributed to the fact that the analytical model fails 
to account for such phenomena. 

The general departure of the theoretical curves from the data led to 


a re-evaluation of the analytical model. It is felt that the weakest 
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Figure 8. Throttling for Normal Injection 
(Unmodified Theory). 
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Figure 9. Throttling For Various Angles of 
Injection (Unmodified Theory). 
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Figure 10. Effect of Mainstream on Jet at High Angles 
of Injection. 


part of the analysis is in the determination of the pressure acting on 
the windward face of the jet by means of the assumption that the press- 
ure varies fromthe static pressure at station 1 down to P at the point 
Of maximum penetration. Although P is fixed since the point of maximum 
jet penetration is the point of minimum flow area for the primary fluid, 
it can be argued that the pressure in the recirculatory region at the 
base of the windward face of the jet may be greater than the Une ansem 
Static pressure. A pressure weighting factor, B, was sought which would 
allow for the mainstream total pressure to act of the windward face of 


the jet at station 1. The following form of 8B: 


po p° S 
=m | fm. |(z) . (14) 
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when substituted into Equation (6) satisfies the boundary conditions: 


BP 


it 
mw 
c+ 
< 
i 
cn) 


gBP = P at ye=h 
The effect of 8 is thus to weight the pressure distribution in such a 
way as to yield higher pressures at the base of the windward face of the 
jet while maintaining the condition P = Pp at y =h. The effective 
average pressure acting on the windward boundary of the jet is now 


expressed with BP replacing P: 


a t f BP dy (6a) 
“| 


If, following previous developments, Equation (6a) is expressed in terms 
of Mach number and substituted into Equation (5), an equation in terms 


of the mainstream Mach number, ih ae appears as follows: 
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Before Equation (10a) can be solved, a suitable value of s must be 
estimated. Figure (11) shows the dependency of g ons. It seems most 
probable that if a stagnation region exists at station 1, the mainstream 
Stagnation pressure would act on the windward face of the jet over a 
short height (y very small) and decrease rapidly to the values calculated 
in the unmodified analysis. From this reasoning, a value of s<<1 would 
be the logical choice based on Fig.11. The modified theoretical curve 
for a = 0° and s = 0.25 is shown in Fig.12. Also shown in this figure 

is the unmodified theoretical curve and the theoretical curve from Nunn's 
| analysis [Ref. 6]. The agreement between the data and the modified 
theory is extremely good up to d,P5/BP = 0.5. For higher values of 
this parameter, the theory predicts greater throttling of the mainstream 


than was found experimentally. 


BP 





Figure 11. Effect of Exponent s on Pressure Weighting Factor,8. 
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Figure 12. Comparison of Theoretical Curves for 
Normal Injection. 
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It should be noted that Nunn's analysis predicts, for a given 
secondary flow, an improvement in the throttling of the primary flow as 
the slot width is reduced. The theory presented in this report does 
not predict such a dependence on slot width. Within the scatter of the 
data presented in Fig 12, it is felt that a dependency upon slot width, 
if any, is not clearly discernable. Further, it should be pointed out 
that Nunn's data included slot widths that were three times greater than 
the maximum slot width examined in this study. 

Figure 13 shows a sequence of oi] flow photographs taken for various 


values of: 


with d, = 0.059 and a = O°. All three photographs exhibit the "Mach- 
bottle" characteristic of underexpanded jets as 1s expected. The photo- 
graphs show very good definition of the windward boundary of the jet up 
to the point of maximum jet penetration. This lends a considerable 
amount of confidence to the assumption that there is a negligible amount 
of mixing between the jet and mainstream prior to the point of maximum 
penetration. 

The shock pattern within the jet plume is not readily discernable 
from the photographs. It can probably be stated, however, that the Mach 
disk occurs somewhere close to the point of maximum penetration. It is 
clearly evident from the photographs that the jet reattaches itself to 
the wall after the point of maximum penetration. At the reattachment 
point the mark of an oblique shock can be seen coming off from the lower 
Surface. This shock serves to curve the jet fluid back so that it 


again flows parallel to the lower block of the test section. 
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c = 0.80 


Figure 13. 0i1 Flow Photographs, d, = 0.059, a Oh 
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Examination of the oil streaks on the surface of the lower block of 
the test section indicated that a small amount of the mainstream fluid 
flowed around the sides of the jet sheet thereby creating a three- 
dimensional problem rather than the hoped for two-dimensional Situation. 

The oil streak which appears in the interior of the jet plume is 
‘thought to be oil which is introduced into the jet at the sides by the 


recirculatory flow in the bubble region. 
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V. CONCLUSIONS 


1. The important parameters in the process of throttling a sonic flow 
by the injection of a transverse jet, are the jet slot width (d,), the 
angle of injection (a), and the ratio of jet to mainstream total 

PRES SUOS . 

2. The sum Mh, at mi, is the controlling factor in determining the effect- 
iveness of the throttling process and not Mm, alone, since the secondary 
fluid, introduced for the purpose of Sipoetano. contributes to the 
overall flow through the nozzle. A practical limit to the throttling 
process is around 60 percent reduction of the mainstream flow. Reduction 
of the mainstream beyond this level requires an increase in the secondary 
flow which is greater than the additional throttling achieved in the 
mainstream. For example, for normal injection, it requires a jet mass 
flow rate of approximately 40 percent of the mainstream reference flow 
to throttle the mainstream by 60 percent. The best overall] throttling 
of a nozzle flow that can be attained is therefore approximately a 20 
percent reduction when the sum of m and m, is considered. Figure 14 
Shows the net throttling attainable (m., 1s m.) based upon the modified 
analytical model. 

3. Within the range of parameters investigated here, the analytical 
model which includes the pressure weighting factor, 8, with a value of 

S = 0.25, is adequate to predict the mainstream throttling up to the 

60 percent level. 

4. To achieve the maximum reduction in mainstream flow for a given level 


of secondary flow, the angle of injection should be a maximum. However, 
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Figure 14. Net Throttling Available for Various Angles 
of Injection. 
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an increase in the angle of injection beyond 30° does not appear to 
appreciably increase the throttling effectiveness. 

5. If jet throttling were used in a proportional amplifier/moment 
producer, the savings in gas expenditure during the no-control phase, 
would be offset by the increase in gas used during the control phase since 
one of the control jets would be on at all times. It is concluded that 
the use of aerodynamic throttling in fluidic control elements is not 
advantageous unless sufficiently higher ol own can be achieved with 
less secondary flow. Suggestions concerning’ future work towards this 


goal are included in the following section. 
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VI. RECOMMENDATIONS 


1. Further study of the throttling process is warrented with regard to 
the effect of the jet slot width. Nunn's analysis [Ref. 8] predicts an 
increase in throttling with a decrease in the slot width, while the 
analysis presented here shows no such dependence on slot width. The 
data in both studies seemed to substantiate the theory presented for 
the analytical model in question. 
2. The analytical model presented here will accurately predict the 
throttling for values of: 

>. Os: 

My 
Further work should be done both experimentally and analytically to 
better understand the throttling process at values of: 

wi 0.5 

My 
3. Further work should be done with regard to the effect of high angles 
of injection on the throttling process. First, an analytical model is 
needed which would account for the fact that the effective angle of 
injection, as seen by the primary flow, is often less than the designed 
for angle of injection. Secondly, an experimental setup is needed that 
would allow for angles of injection approaching 90° with some confidence 
that the angle of injection, as designed into the experimental setup, 
will be felt by the primary flow. | 
4. Regarding the experimental apparatus used for this study, it is 


recommended that if it is desired to use this equipment for future 
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throttling studies, major modifications be made. It 1s recommended that 
the wind tunnel be discarded and that an entirely new test section be 
designed that would mount directly to the transition section from the 

5" pipe to the present test section. In this way, the need for reducing 
the width of the ‘eae by use of inserts could be eliminated and a 
“noticeable improvement in the discharge coefficient should be realized. 
It should be mentioned again that the differential pressure measurement 
required for determining the mass flow rate when using an orifice meter 
was a major source of error. If orifice meters are used in future work, 
it is recommended that they be installed such that the orifice plates 
could be easily changed so that a reasonable pressure drop would be 
measured for all levels of throttling. It is further recommended that 
alternative methods of determining the mass flow rates of the primary 
and secondary systems be sought. Finally, a high pressure compressor 
and air receivers should be obtained to provide for the secondary flow 
system. Although the present system has operated quite well, the high 


expenditure of Nitrogen has made data gathering an expensive proposition. 
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APPENDIX A 


The method of characteristics (MOC) is a well known numerical method 
for obtaining the solution of the complete differential equations of 
motion for supersonic, two-dimensional flow. There are two basic methods 
which use the MOC in analyzing a particular flow field; the Lattice method, 
and the Field method. The field method was used in this study, however, 
a brief description of both methods is included for comparative purposes. 
Both methods assume that the ratio of specific heats 1S constant and 


that the flow is isentropic, inviscid, and irrotational. 


A. LATTICE METHOD 
In the lattice method, the calculation of flow properties progresses 
from a known data line along right and left-running characteristics. 


The flow properties v and 6 behave in a special way along characteristics: 


v + 6 = constant right-running 


v - 6 = constant left-running 

If v and 6 are known along some data curve, then the values of Vp 
and 8p at some point P can be determined from Equations (A-1) and (A-2). 
(See Fig. A-1-a.) 


5 (Sy Wels z (8) - 8) (A-1) 


< 
| 


api : ut 
Once the values for vw and 6 at point P are found, the other flow 
parameters are easily found. It still remains, however, to locate 


point P. 
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(a) (b) 


Figure A-1. Lattice Method. 


The characteristics at point 1 make an angle yp, with the flow 
direction; but if the Mach number changes as the flow progresses down- 
Stream, the angle »p changes and the characteristic lines become curved. 
The curved characteristics are approximated by straight lines or a 
series of straight lines as shown in Figure A-I-b. As the mesh size 
is reduced and more intermediate points are used in locating point P, 


the accuracy of the problem is increased. 


B. FIELD METHOD 

At any boundary change or other flow disturbance, waves are formed, 
and if these waves are weak, they may be used in a manner similar to 
characteristics. It must be remembered, however, that when using weak 
waves as characteristics, the flow properties are discontinuous across 


the wave. 
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In a region bounded by waves, the flow properties are constant. As 
a wave is crossed and a new region (cell) is entered, the flow properties 
change to new values but remain constant in the new cell. The waves 
divide the total flow field into cells of uniform flow. Thus, the 
calculation of properties in the flow field progresses from a reference 
cell to adjoining cells. 

Figure A-2 shows how the jet plume for an underexpanded sonic nozzle 


would be divided into cells by the expansion and compression waves. 







REFERENCE 
CEEL 


Figure A-2. Field Method. 


Expansion Waves - a,b,c,d,e,f,g 

Compression Waves - h,i,j 

Roman numerals - Cell numbers 

All waves are assumed to be of equal strength such that they all 

deflect the flow the same amount &. Expansion waves from the lower 
boundary (a,b, and c) turn the flow away from the center-line by amount 
6 and the expansion waves from the upper boundary (d,e,f, and g) as well 
as the compression waves from the lower boundary (h,i, and j) turn the 


flow back toward the center-line by amount 6. Rotty [Ref. 19] suggests 


that, the cells be identified by four numbers & 4 ; m indicates the 


ay 





number of expansion waves generated at the upper boundary that must be 
crossed in reaching the cell, n indicates the number of expansion waves 
originating at the lower boundary which must be crossed, o the number of 
compression waves from the upper boundary, and p the number of compress- 
ion waves from the lower boundary which must be crossed to reach the 
cell in question. Table A-I shows how the cells in Fig A-2 would be 
identified: 


Table A-I. Cell Identification. 


m oO 
CELL n 


REF. b 4 
. a 
I] (° 0| 
III ( 5) 
IV (3 a 
3 8) 
VI ¢ °) 
VII t 
VIII (3 + 
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If the values for Vp and Op are known in the reference cell, then 
the values for v and 6 can be found in any other cell from the four 


identifying numbers m,n,o, and p. 


v = vp + é(m+n-o - p) (A-4) 


cD 
tt 


Op + é(m - n- 0 + p) (A-5) 


After the flow properties v and 6 are known throughout the flow 
field, it becomes necessary to determine the locations of the wave inter- 
sections so that the flow field can be accurately mapped. To locate the 
wave intersections it must be remembered that the waves are assumed to 
be characteristics and, as such, they make an angle yu with the local flow 
direction 6. The problem now becomes one of finding the directions of 
the right and left-running characteristics for each cell in the flow 
field. The Mach number in the cell can be found from the Prandtl-Meyer 


expression: 


Ge vl TANT! veo (mM? -1) - TANT! Wome =] (A-6) 


The Mach angle, u, is found from: 


w= SIN! @ (A-7) 

and the directions of the characteristics for a given cell are: 
Q@- 4p right-running (A-8) 
at+y left-running (A-9) 


The directions of the waves in Fig. A-2 are calculated in the follow- 
ing manner. The direction of wave (a) is assumed to be an average of 


the left-running characteristics in the reference cell and cell I. 
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Similarly the direction of fave (d) is found from averaging the values 
of the right-running characteristics in cells [I and III. In this manner, 
the directions of all waves can be calculated. Locating point 2 in 
Fig. A-2 is now a simple matter since point | and the vertical distance 
from point 1 to the center-line are prescribed, and the direction of 
wave (a) is known. Similarly, with the location of point 2 known along 
with the directions of waves (b) and (d), the location of point 3 can be 
determined. The calculation of additional points proceeds in an orderly 
fashion from already calculated points along waves whose directions are 
known. 

The accuracy of this method depends on the number of expansion waves 
that are chosen for the initial expansion fan at the nozzle lip. 

The first step in analyzing the flow field for a free jet is to 
determine the initial expansion of the flow at the nozzle lip. The 
flow angle, Av, through which the flow will expand is the difference 
between the Prandtl-Meyer expansion angles corresponding to the nozzle- 
exit Mach number (vy) and the jet-boundary Mach number (vp). Inclusion 
of the flow direction as imparted by the nozzle half-angle Oy completes 
the equation for finding the initial turning angle: 

B= Wy - Up - Oy ~ (A-10) 
This study restricted itself to jets which were sonic at the nozzle exit 
and had a nozzle half-angle equal to zero, therefore Equation (A-10) 
becomes: 
B = -vVp 

The initial turning angle g is equal to the flow direction in cell II 


in Fig. A-2. The boundary Mach number, the only unknown in determining 


= 
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Vp» is a function of the ambient pressure P_, and the jet total pressure 


Pi and can be found from the isentropic pressure relationship: 


ee 
: ~Y 
= (. + hae) (A-12) 


~o | 
c1. O18 


Using the Prandtl-Meyer relationship, as given in Equation (A-6), Vp can be 
determined. 

The initial expansion fan at the nozzle lip, which turns the flow 
through angle 8, is approximated by a finite number of N expansion waves. 


Thus, the strength, or flow deflection, 6, caused by one wave 1s: 


(A-13) 


O 
iN 
=Z|m 


As mentioned earlier, the accuracy of the field method improves as N is 


increased thereby decreasing the effect of any one wave. 


C. APPLICATION OF FIELD METHOD TO FREE JET STUDY 

The field method of the MOC was employed to determine the plume 
boundary of a jet issuing from a two-dimensional slot into a quiescent 
atmosphere. Since this calculation done by hand with any degree of 
accuracy would be a lengthy and time consuming process, a digital 
computer solution was used. 

The computer program presented in this appendix was designed for use 
in the IBM 360/FORTRAN IV system installed at the Naval Postgraduate 
School. 

In the discussion which follows, all procedures are set up for 
calculations below the center-line. Field coordinates are established 
with respect to an origin located on the center-line at the nozzle exit 


ptane. The coordinates are given relative to one-half the width of the 
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nozzle exit; for example, the coordinates of the lower nozzle lip are: 
x = 0.0 and y = -1.0. The computer program maps the flow field down- 
stream of the nozzle to a point where the effect of the first compression 
wave from the upper boundary is felt. Calculations do not proceed beyond 
this compression point. The decision to terminate the program at this 
point was based on the fact that the maximum width of the jet plume will 
occur prior to this point in the flow field. The location of maximum 
plume width is of special interest in this investigation since it has 
been shown by Vick, et.al. [Ref. 12] that the Mach disk occurs close to 
this point. The flow properties undergo a discontinuous change as the 
jet fluid passes through the Mach disk. The MOC will not predict the 
location of the Mach disk, nor will it consider the effect of the Mach 
disk on the downstream fluid. Therefore, further application of the MOC 
downstream of the Mach disk will not give accurate results. 

The only inputs that are required for the program are the pressure 


ratio: 


and N the number of waves chosen for the expansion fan at the nozzle lip. 
The program assumes that the ratio of specific heats is equal to 1.4. 

The first part of the program calculates the boundary Mach number, 
the initial turning angle (FTHETA), and the strength of each wave (DELTA). 
Based on the value of N, the total number of cells in the flow field (NC) 
is calculated along with another field parameter NOC. NOC is the cell 
number for which no compression waves have been encountered in reaching 
the cells numbering from 1 to NOC. Figure A-3 shows the cell numbering 


for a flow field with N = 4. 
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Figure A-3. Cell Numbering. 


Since the program does not treat the flow field beyond the first 
compression wave from the upper boundary, Equations (A-4) and (A-5) 


appear as: 


v = 6(m tn - p) (A-14) 


6 = é6(m- n+p) 


The above expressions are valid only for a jet which is sonic at the 
nozzle exit (vp = 0.0), and with the nozzle half-angle equal to zero 
(6p = 0° ). Equations (A-14) and (A-15) appear in the program as: 

V = DELTA(UE + LE - LC) 


THETA = DELTA(UE - LE + LC) 

where UE is the aonb of expansion waves from the upper boundary which 
are crossed, LE is the number of expansion waves from the lower boundary 
which are crossed, and LC is the number of compression waves from the 


lower boundary which are crossed in reaching the cel] in question. 
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The following six sections of the program assign the appropriate 
values of LE, UE, and LC to each cell, and then calculate v and 6 for 
each cell. After the program has calculated v and 6 Equation (A-6) 1s 
solved iteratively for the Mach number in each cell. With the Mach 
number known, the Mach angle p, and the directions of the left and 
right-running characteristics can be calculated by using Equations (A-7) 
through (A-9). 

At this point the program has stored all necessary information to 
locate the intersections of the expansion and compression waves. The 
intersections of the expansion waves with their own reflections from 
the center-line and the intersections of the expansion waves with the 
center-line are calculated first. Figure A-4 shows the numbering 


scheme for the intersections just mentioned. 


. 1 7 P, 


= Figure A-4. Numbering Intersections of Expansion 
Waves. 
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The location of the wave intersections is carried out in the program in 
the same order as the intersections are numbered. As already mentioned, 
this section of the program locates two different types of intersections: 
expansion waves with the center-line; and expansion waves with their own 
reflections. Each type of intersection requires a slightly different 
procedure for locating the point of intersection. To begin with, ce 
angles, A and B, are defined for the flow field. A is the angle that 
the waves from the lower boundary make with the center-line direction, 
and B is the angle that the waves from the upper boundary make with the 
-center-line direction. Angle A is determined by averaging the directions 
for the left-running characteristics on either side of the wave in 
question. Similarly, angle B is found by averaging the directions of 
the right-running characteristics. Figure A-5 illustrates the method 


for locating an interior point such as point 3. 





Figure A-5. Locating Interior Point. 


60 





Angles A and B are found and then expressed in terms of their tangents: 


a = ERG : LR(II) 


g = RR(I) : RR(V} 


Yn-Y 
TAN A = 4 (A-16) 
ja . 
Ya-¥ 
TANB = 4 (A-17) 
ae 


Solving Equations (A-16) and (A-17) for the unknowns, Xz and y, yields: 


Y5 - y, + X, TAN A - x, TAN B 
X3 = 16 AGE | eee (A-18) 
TAN A - TAN B 


Y3z = (xX, - X5) TANB + yo (A-19) 


Equations (A-18) and (A-19) may be used for locating any interior point 
in the flow field simply by using the appropriate subscripts. With the 
location of point 3 known, the center-line point (6) can be found as 
Shown in Fig. A-6. 

By reine the procedures just described, all center-line points and 
interior points in the expansion region of the flow field can be found. 
The next section of the program locates the intersections of the 
expansion waves with the compression waves and the jet boundary. A new 

numbering scheme is employed for this section of the program. (See 
Fig. A-7.) The previously located intersections along the last expansion 
wave from the nozzle lip are renumbered from 1 to N + 1 and then the 


numbering continues along the compression waves. 


6] 








aX = FAVA 
X6 = X3 + AX 
46 = 0.0 


Figure A-6. Locating Center-Line Point. 


ww 


Figure A-7. Numbering Intersections of 
Expansion/Compression Waves. 
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Locating the interior points in the compression region is done in the 
same way aS in the expansion region, using Equations (A-18) and (A-19). 
The only change that need be made is in determining the value for angle 


A when locating a boundary point, see Fig. A-8. 





Figure A-8. Locating Boundary Point. 


For this case, angle A is the same as the flow direction 6 for the cell 
in question. By calculating angle A in this manner for boundary points, 
Equations (A-18) and (A-19) can be used. 

There are two sections to the output from this program. The first 
is a listing of all the cells in the flow field along with their 
corresponding values of LE, UE, LC, 6, v, Mach number, uy, and directions 
of right and left-running characteristics (RR) and (LR). The other 
section of the output is a listing of the calcualted points along the 
boundary of the jet plume. The boundary points are nondimensional ized 
with respect to the slot width of the jet. 

As was mentioned previously, the MOC will not give accurate results 
downstream of the Mach disk. Since it has been shown that the Mach disk 


occurs close to the point of maximum plume width, it is felt that the 
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boundary points, as located yeti program, are accurate only up to the 
point of maximum plume width. It must be remembered that the field 
method neglects the existence of the intercepting shock in the jet plume 
prior to the Mach disk. It is assumed that the expansion waves from the 
upper boundary pass through the intercepting shock undeflected. Although 
this will result in some error, it has been shown by Vick, et.al. [Ref. 


2) ubhat for: 


20% 


1A 


P. 
a2 
fe 
very close agreement exists between the experimentally determined jet 
boundaries and those boundaries predicted by the MOC after ignoring the 
effect of the intercepting shock. Therefore, it is felt that the program 


will accurately locate the jet boundary up to the point of maximum 


plume width for jets which have a pressure ratio: 
P. 
Sul 
ime 


no greater than 20. 
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